﻿* Encoding: UTF-8.

RECODE PRE_5_combined_model (Lowest thru 0.0001=0.0001) (0.9999 thru Highest=0.9999) (ELSE=Copy) 
    INTO p_trim.
VARIABLE LABELS  p_trim 'p_trim'.
EXECUTE.

COMPUTE logit_p=LN(p_trim / (1 - p_trim)).
EXECUTE.

BOOTSTRAP
  /SAMPLING METHOD=SIMPLE
  /VARIABLES TARGET=MCID_Achieved INPUT=logit_p  
  /CRITERIA CILEVEL=95 CITYPE=PERCENTILE  NSAMPLES=2000
  /MISSING USERMISSING=EXCLUDE.
LOGISTIC REGRESSION VARIABLES MCID_Achieved
  /METHOD=ENTER logit_p 
  /PRINT=GOODFIT CI(95)
  /CRITERIA=PIN(0.05) POUT(0.10) ITERATE(20) CUT(0.5).

COMPUTE Bier=(MCID_Achieved - PRE_5_combined_model) ** 2.
EXECUTE.

BOOTSTRAP
  /SAMPLING METHOD=SIMPLE
  /VARIABLES INPUT=Bier 
  /CRITERIA CILEVEL=95 CITYPE=PERCENTILE  NSAMPLES=1000
  /MISSING USERMISSING=EXCLUDE.
DESCRIPTIVES VARIABLES=Bier
  /STATISTICS=MEAN STDDEV MIN MAX.
